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Abstract 

We derive the universal collapse law of degree 1 equivariant wave 
maps (solutions of the sigma-model) from the 2+1 Minkowski space- 
time, to the 2-sphere. To this end we introduce a nonlinear transfor- 
mation from original variables to blowup ones. Our formal derivations 
are confirmed by numerical simulations. 



1 Introduction 

In this paper we investigate the phenomenon of collapse of degree 1 equivari- 
ant wave maps (solutions of the a-model) from the 2 + 1 Minkowski space- 
time, ]R^+^, to the 2-sphere, S"^. Besides of purely mathematical interest, the 
study of the blowup phenomena for such maps is motivated by the recent 
efforts to understand the singularity formation in general relativity [30] . 

A wave map, <f, from a (rf + l)— dimensional space-time, M, with a metric 
1] to a Riemannian manifold, N, with a metric g, is a critical point of the 
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action functional, which, in local coordinates, has the form 

(This action functional is also known as the a— model). Critical points of 
S{{p) satisfy the Euler-lagrange equation 

where Tjjfj{ip) is the Christoffel symbols on A^. This system of nonlinear 
PDEs is Hamiltonian, and in particular has conserved energy, S{(p), and 
scale invariant in the sense that if (p{x) is a solution then so is ip{Xx). The 
energy, S{(p), is transformed under scaling as 

where f\{x) = ip{\x). Thus the case = 2 of interest for us is the energy 
critical. 

For a map, (p, to have finite energy, it should converge to a constant at 
infinity. In this case for each moment of time, t, if can be extended to a 
continuous map from S'^ to taking the point at infinity to the limit of 
(p{x) at the spatial infinity. Then one can define the degree, degip, as the 
homotopy class of as a map from S'^ to A^. This degree is conserved under 
the dynamics generated by the Euler-Lagrange equations above. 

In our case, M is the 2 + 1 Minkowski space-time, M^+-^ and is the 
2-sphere, with the standard metric g := du^ + sin^ ud9^. In this case the 
degree of (p is an integer (the degree for maps from S'^ to S'^). Moreover, one 
has the Bogomolnyi inequality ([^) 

^if) > 87r|deg(/?|, 

which leads to self-dual and anti-dual equations for the minimizers of the 
static energy for fixed degrees. These equations have explicit solutions (har- 
monic or anti-harmonic maps). These solutions will be written out below. 

Among the maps of the degree k the simplest, most symmetric maps 
are the 'radially symmetric' or equivariant maps which are of the form 
ipk{p,(p,t) = {uk{p,t),k(f)), where {p,4>) are the polar coordinates in and 
[if, 6) are the spherical coordinates in S"^. Then the Euler-Lagrange equation 
for if reduces to the equation: 



u = i\u — — r sin(2w) 

9 n-^ 



(1) 



WM, September 15, 2009 



3 



for Uk- Here A is the 2D spherical Laplacian. Moreover, deg^^fc = Q{uk) = k, 
where 

Q{u) := -{u{oo)-u{0)). 

TT 

Numerical studies of Eqn ([T]) led to a conjecture that large-energy, degree 
one initial data develop singularities in finite time and the singularity forma- 
tion has the universal form of adiabatic shrinking of the degree-one harmonic 
map from to S"^ [4j. Later, it was shown by Struwe [3^ that the existence 
of a nontrivial harmonic map is in fact the necessary condition for blowup 
for 2 + 1 equivariant wave maps. In this paper we address the question of the 
dynamics of the blow-up process. We show that there is < < oo such 
that, as t — > t*, we have on bounded domains in 

^p,t)^f/(p/A(t)), 

where U{p) is the profile of the degree 1 equivariant, static (and in partic- 
ular harmonic) map, minimizing static energy (see below), and the scaling 
parameter X{t), satisfies the following second order ODE: 

XX = T-j^, with a = (1.04)2e~2 ^ .146. (2) 

We expect that, proceeding as in [5], we can show that the error term in the 
above relation is O(A^). 

Note that Eqn shows that if X\t=o < 0, then A < 0, A > for t > 
and therefore A and |A| decrease as t — >■ t,,. Since A^ is the small parameter 
in our analysis (adiabatic regime), our approximation improves as t — *• t^,. 

An approximate solution of Eqn ([2]) with two free parameters (constants 
of integration), t* and c, is (see Section [H] below) 



Mt* -t) = Ae^°'''(x) + c^e'/' 



l_$(_l/2 + lni/2(^)) 
A 



(3) 



where $(a;) = erf(a;) is the Fresnel integral [15]. An exact solution of Eqn 
([2]) is obtained in Section E] (see Eqn (17T]) ). A comparison of the leading 
term of this solution with a numerical solution of Eqn ([T]) is given in Fig 
1. This figure shows that the two resulting curves are indistinguishable for 
times sufficiently close to the blow-up time. 
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° numerics 
f(x) 
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Figure 1 : For a numerical solution that blows up at time t* we plot y = ln ^ ^ ^ as 
a function of x = — ln(t* — t) (circles) and compare it with the analytic formula y = 
f{x) = — ln{a) — \/x + 6, where a = 0.146 and 5 is a free (non-universal) parameter. 
Fitting b we get an excellent agreement between numerical and analytical results. 



Observe that like Eqn ([T]), Eqn ([2]) is a Hamiltonian equation. Its La- 
grangian is 

L := h{X) -\nX, (4) 

where the function h is defined by h"{x) = —l/f~^{x^/a) with f{x) = 
a;ln(l/x) (see Section [7j). 

The local well-posedness for the wave map equations in Sobolev spaces 
was proven in [T71 [HI [19], while the global well-posedness for small initial 
conditions, in [20l [36l ill iSl iQl E] (see also [T0l[IIl[T2l[2ll[22l[33l[3ll[36l 
dU HSj). The research on the problem of blowup for the wave maps started 
with numerical work [H [2S1 • (We do not review here related works for 
nonlinear wave equations.) 
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The first numerical evidence for singularity formation for 2 + 1 equivariant 
wave maps to the 2-sphere was given in [1]. In this paper (concerned only with 
k = 1 homotopy) the authors showed that blowup has the form of adiabatic 
shrinking of the harmonic map and formulated conjectures about blowup for 
large energy, blowup profile and energy concentration and that A [t) /{T — t) 
must go to zero. As was already mentioned, it was shown rigorously in 
[27| that the existence of a stationary solution is a necessary condition for 
the blowup to take place. The blowup scenarios were further numerically 
investigated in fi6[ [26] (see references therein for additional works) . 

The first rigorous results on the blowup rate and profile were obtained in 
[211 [21] • In particular, [3T] has obtained the lower bound on the contraction 
rate X{t) for k > 4 wave maps. As it turned out this lower bound conforms 
exactly to the dynamical law derived for the 4 + 1 Yang-Mills k = 1 equiv- 
ariant solutions in [5j using a formal but careful analysis, explained below 
in this introduction, justified by numerical computations. (Earlier numeri- 
cal analysis for the latter model was announced in [6] and described more 
completely in the survey [2].) (It was noticed in [3l| (see below), that the 
k >2 wave map equation are similar to the 4 + 1 Yang- Mills one for k = 1. ) 
Finally, for each b > 1/2, [2l| has constructed special solutions of the k = 1 
equivariant wave map equation, Eqn ([1]) with k = 1, which blow up at the 
rate A(t) ~ (T-^)l+^ 

Equation ([T]) belongs to a general class of semilinear wave equations in 
M^+i of the form 

u = Au~^fiu). (5) 
P 

In the case of f{u) = ^sm{2u), Eqn ([5]) is, as was already mentioned, 
the equation for the profile of the equivariant wave map from the 2 + 1 
Minkowski space-time of degree k, R^+^, to the 2-sphere, 5*^. More generally, 
dSD is satisfied by equivariant maps for the case when N is the surface of 
revolution with the metric g := du"^ + g'^{u)d6'^, where g{u) is related to f{u) 
as fiu) = g{u)g'{u). 

In the case of f{u) = 2u{u'^ — 1) the corresponding equation, 

2 

ii = Am - — (m^ - i)m, (6) 
P 

is related to the equation for equivariant Yang-Mills fields of degree 1 in the 
4 + 1 dimensions. 
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Note that 

(i) Eqn ([5]) is invariant with respect to the scahng transformation, 



A' A 



(ii) Eqn ([5]) can be presented as a Hamiltonian system with the standard 
symplectic form and the Hamiltonian 



with F'{u) = f{u). The scahng properties of the Hamiltonian H{u,v) imply 
that the dimension d = 2 is the critical dimension for Eqn ([5]). This is the 
dimension treated in this paper. 

We assume now that f{u) is a derivative of a double-well potential F{u), 
i.e. F{u) is nonnegative and has at least two global minima, say at a and 
b for some b > a, with F{a) = F{b) = 0, and no minima between a and 
b {F{u) = I sin^ u and a = 0, 6 = vr in the case of f{u) = Ysin(2M) and 
F{u) = — 1)^ and a = —1, 6 = 1 in the case of f{u) = 2u{v? — 1)). In 
this case Eqn (jS]) has the following features: 

(A) For each A; G N, Eqn ([5]) has static solutions, f/fc(p) and U^k{p) = 
—Uk{p); they have topological degrees Q{Uk) = k and Q{U-k) = —k; 

(B) For k = 1 , the solution Ui{p) is monotonically increasing from a to b, 
while U-i{p) is monotonically decreasing from a to b; 

(C) The solution, Uk{p), is a minimizer of the static energy functional E{u) 
under the constrain, Q{u) = k, on the topological charge; 

(D) Eqn ([5]) conserves the topological charge Q{u) := ^(u(cx3) — u(0)). 

Existence of the solutions Uk{p) follows from the Bogomolnyi argument, 
see above. The solutions Ui{p) and f/_i(p) are called the kink solution and 
anti-kink solution, or simply kink and anti-kink, respectively. Since the anal- 
ysis for k < can be obtained from analysis for the case > by simply 
flipping the signs, in what follows we assume that k > 0. Note that though 
Eqn is scale invariant, its static kink solution Uk{p) are not. Hence Eqn 

has an entire family, Ukij), of kink solutions (symmetry breaking). 





WM, September 15, 2009 



7 



From now on we concentrate on the kink solution, Ui{p) and omit the 
subindex 1: Ui{p) = U{p). 

There is a feature of Eqn ([5]) which is not apparent at the first sight but 
which plays an important role in our analysis of the collapse. The fact that 
the kink, U{p), breaks scale invariance manifests itself in appearance of the 
dilation zero mode 

C(p) := lpd,U{p). 
This is a zero eigenfunction, Lp( = 0, for the linearization operator 

(negative Frechet derivative) for the r.h.s. of ([5]) around U{p). This zero 
mode presents an obstruction to solving Eqn (2) perturbatively, starting 
with U{p), which can be resolved by a modulation theory, provided ( is an 

function, i.e. one can use a Hilbert space spectral theory. 

Thus equations of the form ([5]) can be organized in two classes according 
to which of the following two properties takes place 

(i) ( is in 

(ii) ( is not in L^. 

The Yang-Mills equation, ([6]), belongs to the first class while the wave map 
equation, ([T]), belong to the second. Indeed, the kink solution for ([6]) is 
U{p) = and the corresponding zero mode is ({p) = ^^^^^2)2 (see [21 E])- 
For Eqn ([T]) with k = 1 the kink solution is 

U{p) = 2 arctan p, 

while the scaling zero mode is 

C(p) := ]-pdMp) = 

2 1 + 

Clearly, ( is in the former case and is not in the latter case. (This 
is possible due to the fact that the operator L has no spectral gap: cr(L) = 
acont{L) = [0, oo) (see below). The fact that there is a problem with the mod- 
ulation approach due to the nonintegrability of the zero mode was pointed 
out by P.Bizoh in 2001, [BJ.) 
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Note that Uk{p) = 2arctan(p^) and the corresponding zero mode is 
square- integrable for k > 1. Thus in this case we expect that at least the 
formal analysis of [5] of the Yang-Mills equation should go through. (Higher 
degree equivariant, static solutions for the Yang-Mills equations in 4 -|- 1 
dimensions, known as instantons, can be found in [ll [32].) 

We are interested in solution with initial conditions near the kink manifold 

Mk^nk■.= {Uip/X)\\>0}. 

The fact that the square integrability of of the zero mode ( plays an important 
role in analysis of such solutions can be gleaned from the observation that the 
effective action S{Ux) on the family Ux{p) := U{p/\) (the 'effective action' 
of A) is equal to 

siux) = l Jmncr-v{u)}dt, 

where V{u) := ||VmP + ^F{u), with F'{u) = f{u), and diverges, if C ^ 
(For a connection to the geodesic hypothesis see [STJHS]). Here and in what 
follows A = §. 

We present heuristic arguments motivating our approach. It is natural to 
guess that for an initial condition close to the manifold Mkmk the solution 
evolves along this manifold. Let U{^^) be the projection of the solution on 
this manifold. If for this projection A(t) ^ as t ^ for some t*, then the 
solution collapses at the time t^,. With this in mind we look for solutions to 
Eqn ^ in the form 

u{p,t) = U{x) +w{x,t), (9) 

where x = p/A, a blow-up variable, with A a slowly varying function of time t 
(we do not pass to the blow-up time variable). Note that while in a standard 
approach the scaling. A, is fixed at the very beginning (with corrections 
at certain scales possibly considered later on) we leave it free and we look 
for a differential equation for A which guarantees that \w\ <^ 1. However, 
this simple procedure which works in the case of the Yang-Mills equation 
mentioned above (see [3 [31]) does not work in the present case as we explain 
below. 

Note that if G L^, then A(t) is uniquely determined by the orthogonality 
condition 

(C,w) = o. (10) 
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U ( ^ L'^, then this condition is not well defined unless we assume w belongs 
to a space of sufficiently fast decaying functions. 

Substituting decomposition IQ into ([5]) leads to the equation for the func- 
tion w and parameter A: 

L^w + F{w) = -X'^, (11) 

where F(w) absorbs higher order terms {F(w) = A^^|^ + N{w), N{w) = 
nonlinearity in w) and is the linearization operator for the r.h.s. of ([T]) 
around U{x) given by ([H]). The operator is self-adjoint. The scaling zero 
mode, (, is a zero mode of this operator: L^C = 0. Since ( is positive and 
not we conclude by the Perron- Frobenius theory that a^L^) = [0, oo) and 
is not an eigenvalue of L^. 
We compute explicitly 

A'^ = A2[-29i(AA-i)C + 2{XX~'fxd,C] 

= -2AAC + 2A2(C + xa^C)- (12) 

We multiply Eqn ( ITTi) scalarly by C{x). Though ( is not one can show using 
a limiting procedure that {(, L^w) = 0, provided w = o{x) and d^w = o(l) 
at oo. Thus we obtain 

A^(C,^) + (C,W) = 0. (13) 

Following [5j we try to develop a perturbation theory in the small pa- 
rameter A^ assuming that term AA is of the order o(A^) (and A < 0) and 
similarly for higher order time derivatives of A, e.g. dt{XX) = O(A^), etc. 
Furthermore, if our assumption that \w\ <^ 1 is correct and the integral in 
(C, F{w)) is convergent, then we can drop the term {(, F{w)) in f[T^ . Hence 
we obtain in the leading order O(A^) 

X'{C,C + xd,O = 0- (14) 

Considering the integral on the l.h.s. over a bounded domain and integrating 
by parts one shows that the inner product on the l.h.s. is 

1/2 lim(a;^C^(a;)). (15) 
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For Eqn ([6]) this is so we can solve Eqn (fTTl) in the leading order, w = 
—X^L^^{C + xdxQ- Plugging this result into Eqn ( fT3ll and keeping only the 
terms up to the order 0{\^), we obtain the equation for scaling dynamics, 

AA = ^A^ (16) 

in the leading order O(A^) (see [5]|3T]). Next, in order to obtain a correction 
to this equation, we use f|T3|) at the order 0(A®) to solve Eqn ffTTl) to the order 
O(A^) and plug the result to f|T3l) . However, at this step we run into loga- 
rithmically divergent terms. To overcome this problem we use a multiscale 
expansion, by introducing an additional scale at infinity (see [S]). 

For Eqn ([1]) with = 1 we have lim^^oo(3^^C^(3^)) = 1 cind so we go to the 
next term, — 2AA||C||^, and discover that it diverges logarithmically. Thus for 
Eqn ([1]) with fc = 1 one runs into a problem right away. This shows that 
decomposition ([9]) is incompatible with the condition \w\ <^ 1. 

The problem for Eqn ([T]) with k = 1 mentioned above can be also seen 
in a different but related way. Let us try to solve Eqn ffTTl) by perturbation 
theory. In the leading order we drop the term F{w) to obtain the leading 
order approximation to the solution: w = L~^Lp, where = — A^^^ and 
is understood as the Green function of the equation Lw = ip (see Section [3]). 
It is easy to check, using Eqns fl55]) - flMl) of Section [2] below, that if A 7^ 
then the function L~^(f grows at infinity as xlnx, and a straightforward 
perturbation theory fails. (Not only the correction w is large at 00, its 
energy is infinite.) 

The point here is that the function U{p/X) is not a good adiabatic solution 
to Eqn with k = 1: 

X\d^U{x) - A,U{x) + sin(2f/(x))) 

Zp 

= -2AACM+(^. (17) 

where x := p/X and where we used ( fT2l) and the relation ( + xdxC = 2x(l + 
^2^-2^ The r.h.s. is not L^. The problematic term is 2XX({x). In particular, 
it leads to the logarithmically divergent term, 2AA||C||^ in the orthogonality 
condition. Hence, one has to find a better leading term. 

We deal with the problem above by introducing instead of the linear, 
one-parameter transformation, p ^ p/A, a nonlinear, three-parameter trans- 
formation, p /(p. A, a, P), chosen so that U{f{p, A, a, P)) becomes a better 
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approximate solution to Eqn ([T]) with k = 1 than U{p/X). In particular, the 
problematic term 2XX( entering the r.h.s. of Eqn (JT71) is canceled and there- 
fore the large p divergence in Eqn f|T3|) mentioned above is eliminated. Thus, 
instead of ([9]), we look for solutions of Eqn ([1]) in the form 

u{p,t) = U{y) + w{y,t). (18) 

We consider initial conditions close to U{y) = U{f{p,X,a,P)) (we do not 
specify the norm, the latter must be determined by a rigorous analysis, see 
e.g. [SI])- After this we proceed as above with Eqn The conditions 

<^ 1 and w — at p — >• cx) and constraints on the energy and its 
fluctuations imply the differential equation ([2]) on the parameter A = X{t). 
We expect that proceeding as in [5J one can obtain corrections to Eqn ([2]). 

The paper is organized as follows. In Section 2 we introduce a change of 
variables, p f{p, A, a, f3), depending nonlinearly on the original variable p 
and on the scaling parameter A~^ (and depending on additional parameters 
a,P). This is our main new idea. In Section 3 we derive, modulo some 
technical details which are provided in Appendices 2 and 3, an approximate 
solution to Eqn ([1]) with k = 1. In Sections 4 and 5 we use an orthogonahty 
condition of the type of ffTOl) . the smallness condition on energy fluctuations 
and the minimum condition on the energy of the approximate solution in 
order to find our main equation on the scaling parameter A, Eqn ([2]). In 
Section 6 we find exact and approximate solutions of Eqn ([2]) and in Section 
7 we show that this equation is Hamiltonian. In Appendices 1-5 we provide 
technical calculations used in the main text and explanations of the numerical 
approaches. 
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2 Nonlinear blow-up variables 

In this section we introduce a nonhnear, three-parameter (scahng) transfor- 
mation, p f{p, A, a, P) of the independent spatial variable p. This replaces 
the standard, linear, one-parameter transformation, p ^ p/X. We denote 
X := p/X and define y = y{x, XX, a, [3) as the solution of the equation 

AA o / \/ XX ^-n\ T ^ 

y = X —X In — —y x , n x < Xcr 



2 V P 

?/ = 2?/er + Ya;Mn|^^^?/V~°^, if x > x^^, (19) 

where < a < 1, /3 > 0, Xcr = Xcr{X,a,P) and ycr = ycr{^,Ci,P) solve the 
equations ^(r.h.s. (fT9|) ) = and the first equation in (fT9!) . We can write 
Xq^ and ycr as 

^cr=(^) ' (3 ln7 + l-a)-i/2 (20) 



and 



where 7 = 7(0, /5) = ^^'-^y'^x^ " is a solution of the equation 



_ 2 21n7 + l - a 

^■^^"'.mJ (31n7 + l-«)3/2' ^^^^ 



7 (31n7 + l-«)^/2+" ^22) 



^ (21n7 + l-a)" 

Eqn fl22l) is well defined for In 7 > — 1(1 — a) and in this range it has a unique 
solution for each a and /3. We denote per '■= Xxcr- 

For X < Xcr the r.h.s. of f|T9|) is decreasing from 00 to —00. Hence the 
equation flTI?]) has a unique solution for x < Xcr- Furthermore, for x > Xcr 
the r.h.s. of f[T^ increases logarithmically in y from — cxd to 00 and increases 
monotonically in x. Since for x = Xcr (fT9|) has a unique solution, it has 
exactly two solutions for x > Xcr- Of these two solutions we choose the 
greater one. 

Finally, we notice that the function y = y{x, XX, a, (3) increases monoton- 
ically in X for X > 0. Indeed, for < a < 1, the equations ^(r.h.s. (fT9l) ) = 
and f|T9|) have a unique solution (x = Xcr, y = ycr), and therefore so is the 
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equation = 0. Since y is obviously increases monotonically in x for x 
sufficiently small and sufficiently large, it does so for all x > 0. 

Write v{y,t) = u{p,t), where y = y{x, XX,a, (3) is given in Eqn f|T9l) . In 
the new variables, Eqn ([1]) with k = 1 becomes 

d'^v 1 dv sin(2f ) ^ , , , . 



dy'^ y dy 2y 



*W - 51 f^X + X^l S . f ^ . 1^1 1^ - A=S'^ (24) 



where 



y^ {\ X J dy"^ \ X dx J dy dt'^ 

and d'^/dt'^ is the total derivative in t (i.e. taking into account that y depends 
on t). Here the function x is defined according to the equation 

1^ = ^ + X. (25) 

dx X 

Eqn (123!) is our transformed equation. 

Initial conditions for (|23ll are chosen to be close, in an appropriate norm, 
to U{y) = U {y{x, XX, a, P)), where, recall, U{p) is the static - kink - solution 
to Eqn ([T]). To simplify the exposition we take the initial condition to be 
just U{y{x,XX,a,l3)). 



3 Approximate solution of Eqn ( 1231 ) 

Let y = y{x, XX, a, P) be the transformation defined in the previous section. 
We look for solutions of Eqn fl23|) in the form 

viy,t) = Uiy)+wiy,t), (26) 

where w is a small correction. We plug this decomposition into Eqn fl23l) to 
obtain 

LyW + N{w) = <iJ{U + w), (27) 

where operator Ly is defined in Eqn (|8]), N{w) is the nonlinear in w term 
defined by this equation and the function \E'(f) is defined in fl2^ . To find an 
approximate solution of the latter equation we drop the nonlinearity, N{w), 
and the term w in "^{U + w) to obtain the leading order equation 



LyW = ^, 



(28) 
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where ip{y,t) := tlj{U{y)). The latter function is given expUcitly by 



1/2 [ x(l + ?/2)2 I -\- y"^ ydx xj (l + ?/2)2 



where the variable x is connected to y through fll9p . The counter-term which 



removes the undesirable term in is dx/dx — 2x1 x. As the result we have 



i,{y,t)=o{}?^^^0^ io.y>y^, (30) 

(see Appendix 3). 

Remark 1. To justify dropping from Eqn 0271) (it enters into ip{U + w)) 
is not a simple matter. In a rigorous approach one looks for solutions of (1271) 
in the form 

v{y,t) = U{y) + W{y)+r^{y,t) 
and shows that the fluctuation rj is small (cf. 131]). 

Recall that the operator Ly entering Eqn fl25]) has a zero mode: 

LyC = 0. (31) 

Multiplying Eqn (1251) scalarly by ( and using the self-adjointness of the op- 
erator L (and some elementary limiting procedure) and (1ST]) , we obtain 



oo 

/ 





dyy = 0. (32) 



This is a (necessary) solvability condition for Eqn (1281) . It gives an equation 
on the parameters A, a and f3. (Note that it is an approximate solvability 
condition for the exact equation (1271) . 

So far we obtained one equation for the three parameters A, a and (3. 
To derive another equation we analyze the approximate solution to (l27j) we 
obtained: U{y) + W{y), where W = L~^ip with tp satisfying (|32|) . Our goal in 
the rest of this section is to isolate the leading contribution to W . This will 
be used in the next section to derive the second equation for the parameters. 
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To find L ^ip we compute the Green's function for the operator L. Two 
hnearly independent solutions of the homogeneous equation Lw = are 



I + y 2, ly I + 



(the first of these solutions is just the scaling zero mode, C, the second solution 
is found in Appendix 1). Hence by the ODE theory 

ry ry 

L'^ijj = cwi + wi W2ipy'dy' — W2 Wiipy'dy' (34) 
Jo Jo 

where c is chosen to guarantee solvability of the equation to the second order 
correction term or by minimizing the energy. 

We find the leading contribution to the solution w = L~^ip of Eqn (l28l) . 
In what follows we use the following assumptions 

< AA < < 1, XdtiXX) = O(A^), Xdt{X^) = O(A^), (3 = 0(1). (35) 

We see from ([20]) that 

Xcr = Pcr/X ~ (AA)-^/2 ^ i_ 

Consider first the region y < y^r- In this region if) is given by ([7]), Ap- 
pendix 3. The latter equation shows that for y -C ycr the leading part of -0 
is 



3 8AA ( In ( + 1/2 - AX^ 

where, recall, x := p/X is connected to y through (fT9l) . Now, let 

^2:=^-V'i- (37) 

According to (1341) . the general solution, W, of Eqn (l28ll decreasing at 
infinity is of the form 



W{y,t) = CiWi{y) + wi{y) / W2is)'ilj{s,t)sds 

Jo 

POO 

+W2{y) / wi{s)ilj{s,t)sds, 



y 
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with wi = ( and W2 defined in (l33ll . The function W2 is singular at y = 0. 
Hence W is bounded only if the condition (l32l) is satisfied. 

The function V'ld/, t) is localized on the scale ?/ ~ 1, decays as in the 
region 1 ^ y <^ y^r and decays as ?/~^ in the region y ^ y^r (though we 
are considering at the moment only the region 1 <^ y <^ y^r, the latter fact 
allows us to extend the integrals to the entire real axis with a small error). 
The function iIj2 is localized at the large scale, y ^ ycr ^ ^- After some 
lengthy computations we find for y <^ ycr 



W{y,t) = ^ I w^{s)^^{s,t)sds + 0{\^\\y^). (38) 



oo 







Using the expression fl36|) . it is easy to show that 

Wi{s)ipi{s,t)sds = O(A^). 



4 Energy of the approximate solution and 
the equation on A 

We compute the energy of our approximate solution u{p,t) = U\^a,p{,p) + 
Wiy), where Ux^a,i3{p) '■= U{y), with y = y{x, XX,a, P) and U defined in 
the introduction, and where W = L'^ip, the solution to Eqn (1251) (see the 
previous section). Due to ([7]), the energy functional is 

E{u) := {-e + -|Vm|2 ^—^^r?u)pdp. (39) 

Inserting the approximate solution into this expression, we obtain that E{Ux^a,f3+ 
W) = E{Ux,a,f3) + 5Ei with 

6Ei:=0{XY,,)(^j^ wi{s)Ms,t)sds^ +0{X^). (40) 

Furthermore, we have that 

E{Ux,a,p) = E{U) + 5E^ with = 0{}? ln(l/A2)). (41) 

We require that the energy correction due to the fluctuation, W , be much 
smaller than the one due to the modulation: 



(42) 
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Since wi{s)ilJi{s,t)sds = O(A^) and i/cr = 0(— ^) ^ 0{-i-), this implies 



that the integral in the leading term in the above expression for SEi must 
vanish: 



/ 



dyyCiIJi = 0. 



(43) 



This gives an implicit equation on the parameters A, a, (5. 

In the leading order, we can replace y hj x = p/\ (see the first equation 
in (JT9l)). so that Eqn ( H3i) becomes 



oo 

/ 



dxx- 



X 



iXXx 



1 + X^ [ [1 + X 



2^2 



In 



AA 



P 



-X 



1/2 



+ 



AX^x 

(1 + X2)2 



(44) 



Computing the integrals in (jUj) (see Appendix 2 for detailed computations), 
we obtain 



A^ + 2AA 



.nl-^,.l 



0. 



(45) 



This is our explicit equation for the parameter A. It depends on the additional 
parameter (3 whose value we still have to determine. Since in the leading 
approximation [y x) the first equation on the r.h.s. of flTI?]) is independent 
of a, then so are the resulting equations flUl) and fH51) . Eqns (jlSD and ([2]) 
coincide, provided 

a = P^ e-l (46) 

Clearly, solutions of Eqn (H5ll have the property (135|) assumed above. 
Moreover, if A(0) > 0, A(0) < 0, then, by Eqn (jlsD, A(t) > 0, X{t) < 
0, X{t) > and therefore A(t)^ < A(0)^ for t > 0. As t ^ t*, |A| decreases 
so that our approximation improves as t — > t,,. 

Thus it remains to find the value of the parameter (3 . To this end we use 
the condition fl32l) and minimization of the energy of the leading part of the 
approximate solution. 



5 Values of the parameters a and (3 



In this section we derive an equation on the parameters a and (3 and use this 
equation together with the energy minimization to find the values of these 
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parameters. We assume ( ]35l) and that, at least in the leading approximation, 

a,P are independent of t. (47) 

In what follows we do not display the dependence of the quantities involved 
on XX (however, the dependence on X is displayed). 
First note that Eqns (l32l), ([37]) and (iH]) imply that 

f^Vcr 



ip2dy + / ipdy = 



(48) 



From Eqns (12T!) and fj45l) we obtain easily 



and 



1 ^^(AA) 
AA dt 



1 ^ dXfyp 



2X\n-^ ( ^ 
AA 

-Aln-Mi 
,AA 



1 + 1/ln^ 



1 

AA 



o(i/H 



(49) 



As a result in the main approximation in 1/ ln(l / AA) we should keep in the 
expressions for the functions ip and ^p2 in Eqn fHHl) only terms proportional 
to A^. The latter terms are given in ([7]) and fl96|) in Appendix 3. The 
most important region in the above integral is where y is of order of ycr 
{ycr = 0(— ^) ^ 1 due to Eqn ( !2T1) and the condition (l35l)). As a result we 



AA 



can neglect 1 compared to y^ in (jTl) and ( !96l) . Then Eqn ( l48l) reduces to the 
equation 

I{a,P) ■.= hia,P) + l2{a,P) = 0, 



where 



ycr 



dyx^ 



XXx^ In 



/3 



— a 



2(1 -r^ 



(50) 



2(1 -r) 



AAx / ax^ 



y[} + 

y2 



2j/ 



1/1 + 



lAAxfV 
2y / 



xil + 



2y ' 



Qax 



3(1 



r 



1 _|_ aAAx£^ 



2 1 



V ) 



3 In 



2?/ ) 

13 y 



1 + 



2?; 

+ 1 - a 
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and 



dy 



XXx 



ax 



X 



y2 



2Y 



6ax 



2y ' 
3(1 



a) 



2 1 + 



y 



2?/ ; 



3 In 



aXXx^ 



13 y 



+ 1 - a 



2y ) 



Here 



Y = 1 



XXx^ 



3 In 



/XX 



1 



a 



(52) 



(53) 



One can further evaluate Ji and I2 by changing the variable of integration 
y in (ISTl) and ([5]) to 2; as 



Z = X/ Xcr = pI Pc 



(54) 



where a; and functions of y, AA, a, (3 are given in (1191) and the defini- 

tions following this equation, and compute the resulting integral numerically. 
In particular, one can show that for a = 0, Ji(a; = 0, (3) = 1, 12(0; = 0, (3) = 
and therefore I{a = 0,(3) = 1, independently of the value of (3. Thus we 
cannot take a = in our transformation f|T9|) . 

We chose the parameters a and (3 which minimize the energy E{a,(3) : = 
E{Ux,a,i3), where, recall, Ux^a,/3{p) '■= U{y{x, XX,a, P)), given that the equa- 
tion (l50l) . I{a,P) = 0, holds. To find these minimizers we use Eqns (l39ll and 
U{p) = 2 arctan p to rewrite the energy E{a,j3) as 



E{a,P) 



dpp- 



l + y 



2\2 




+ 



dy_ 
dp 



+ 



(55) 



We find numerically (see Appendix 5 for the analytical part) that the energy 
E{a,j3) is minimized on the curve I{a,(3) = at the point 



/3o = 1.04 and ao = 0.65436. 



(56) 
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This is a special point for the curve I{a,(3) = 0. Our numerics show 
that while the functions a = and /3 = P{a) determined by the equation 
I{a,(3) = are double-valued, their branches originate exactly at this point 
(and form a wedge there). So the equation I {a, j3) = has a unique solution 
only for /3 = /3o or for a = ao and has no solutions for P > Pq or for a < ao- 

Substituting [3 = (3q = 1.04 into Eqn fH5]) . we obtain the following value 
for the parameter a: 

a = 0.146. 
This proves Eqn ([2]) with a = 0.146. 



6 Investigation of Eqn 

In this section we find an approximate solution to Eqn ([2]) (which is, up to 
a redefinition of the parameters, the equation fj45|) ). Iterating this equation, 
we find, in the leading approximation, the following equation 

I ^ (57) 

A Aln(a/A2) 

Solution of the Eqn fl57|) with two free parameters of integration, c > and 
r, is 

A 

^{t*~t) = j rfa; e'-^'^'^i). (58) 







Changing the variable of integration as ^^(^f J ~ ^"^^ reduce Eqn ( 158|) for 

the parameter A to the form Eqn ([3]) given in Introduction. 

Now we derive an exact expression for a general solution of ([2]). We 
introduce the function f{x) := xln(l/x). For < x < this function has 
the inverse, /~^(x). Using this inverse we rewrite Eqn ([2]) as 
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(Note that for x — > 0, / ^(x) = ^ + so in the leading approximation of 
( I59l) gives (l57j) .) Integrating equation (l59l) gives 



InA = F(A) where F{y) = - -^g{z), 



(60) 



with (7(z) := z// ^{z). Using the equation /(/ ^(y)) = y, or, more exphcitly, 
f~^{y) \n{l/ f^^{y)) = y, we find that the function g{z) satisfies the equation 



Differentiating the latter equation, we find 

9{^) 



(61) 



(62) 



Using this equation we integrate 



dz 



This gives 



g{z) = - / dzg\z){g{z) - 1) = --{g{x) - \f + const. (63) 



(64) 



F{y) = -h,g{y' la) -\f\ const, 



which together with Eqn fl60l) yields 

for some constant c. The latter equation can be integrated as follows 

dx 



(65) 



v^(t* - 1) 







g-Ul + 2^: 



1/2- 



(66) 



Next we find the function g ^{x). The definition of the function f{x) implies 
/(e~^') = xe~^, which yields 
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which, in turn, leads to g{xe ^) = x, which finally gives the expression 



g ^(x) = xe ^. 



Now Eqns m&i and (l68ll imply 



^{t* -t)= dx 



i/2+yR|) 

Changing the variable of integration as ln(c/x) = we find 



oo 

^{t* -t) = V2ce^/^ j 



dz 



We obtain for A/ln(c/ A) ^ 1 the approximate expression 

A 



v^(t* - 1) 



v^[ln(c/A)]V4 



,l/2+^ln(c/A) 



4v/ln(c/A) 321n(c/A) 



(68) 



(69) 



(70) 



(71) 



Eqn ([3]) is an approximation for this exact expression, it differs from the latter 
by a slowly varying factor which can be found in the next approximation to 



7 Hamiltonian Formulation 

Eqn ([2]) is a Hamiltonian system. Indeed, it can be obtained from the Lan- 
grangian 

L = h{\)-\n\ (72) 
where the function h is defined by 



WM, September 15, 2009 



23 



with f{x) = x\n{l/x) (see Section 5). Now the generahzed momentum, 
Hamiltonian and energy can be found in the standard way. In particular, the 
energy is given by 

(9/ 

E = -X^ + L = -Xh'{X) + h{X) - In A. (74) 
dX 

This is the energy conservation law. On the other hand, differentiating Eqn 
f l7^ w.r.to t, we obtain the equation of motion ([2]). 



7. Conclusion 



We presented detailed arguments that for an open set of initial conditions 
close to the degree 1 equivariant, static wave map, the solutions of the wave 
map equation (cr-model) collapse in a finite time. Near the collapse point 
the solutions have a universal profile given by the modified the degree 1 
equivariant, static wave map depending on a time-dependent parameter A. 
This parameter describes the rate of compression (scaling) of the collapse 
profile. We derived a second order Hamiltonian dynamical equation for the 
scaling parameter, A. We also found approximate solutions of this equation. 
These solutions are of a rather complex form. They are in an excellent 
agreement with direct numerical simulations of the wave map equation. 



Appendix 1 

To solve the equation Lw = g, we should find first of all two linear 
independent solution of linear equation 

L^w = 0. (75) 

The first solution of this equation is the scaling zero mode ( 



u.^C^^. (76) 



The second solution W2 satisfies the inhomogeneous equation of first order: 

W1W2 — W2w'i = —. (77) 

X 
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The standard solution of this equation is 

2 1 1 

W2 = Wiz; z' = X -\ \- —; 2; = C + — + 21na; - — r. (78) 

X x-^ 2 2x'^ 

Setting C = 0, we obtain 

X 2a; In X 1 , 

To obtain general solution of the equation L^w = g, we rewrite it as a 
first order ODE 



Two linear independent solutions of A1.5 are 

Wi\ { W2 



w: . (80) 



(81) 



By the method of variation of constants we look for a general solution of 
inhomogeneous Eqn flHOl) in the form 



w 



where ci,2 are functions of x. Inserting (182|) into Eqn (IHOj) . we find 

dci dc2 

— = xw2g; = -xwig. (83) 

ox ox 



Appendix 2 

To derive Eqn fj45|) from Eqn fj44|) we should calculate two simple integrals. 
One of them is 



00 

dy If dx X 1 



1 + 2/2)3 2 J {1 + x 



(84) 
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The second integral is {e ^ 0) 



dy y^lny 1 f dx xlnx 

(1 + 1/2)3 = lj (l+x)3 



-- / \nx d 

4 



x + 1 2 (x + 1 



(85) 



- In £ + - 
8 4 



dx / 1 



X \x+l 2 (a; + 1)2 



-\ne + - 
8 8 



dx 



1 



1 1 

+ 



X X + 1 ' (x + 1)2 

Using the values of these two integrals, we obtain Eqn (145|) from Eqn (jHj). 
Appendix 3 

Now we will find an explicit expression for the inhomogeneous term ijj. 
We consider separately two domains {y < ycr} = {p < Per} and {y > ycr} = 
{p > Per}- First, we compute ^ and 

Recall the notation x := p/X and 



Y ■=!- 



XXx^ 



3 In 



/XX 

T 



+ 1 - a 



(86) 



and let 
A :=ln 



AA 



aXXx^ 1 aAAx^ 
X := IH and Z := 1 (87) 



^y 



21/ 



In the domain y < y^r we have 



a|^-Ax| + 0(x3a|(AA 



A^ 



Y Xx\,d 



2X' -XX]x— + —X 



dt 



XX^A + l-a + ^ 



A2 



X 



+ ( x^AA^ ( AA 
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where 



d_ 

di 



\\3A + 1- a + 



A2 



X 



XXXjax^Y^ 6ax 3(1 - a) 2(2Z - 1) (3A + 1 - a) 



X 



X2 



In the domain y > ycr we have 



(89) 



2A^.-AAx)--— A3- 



AA3A+l-«-f 



+0 x^AA— AA 



where 



d_ 

dt 



dt 



A2 Z 



AAA f F2 ^ 6ax 3(1 - «) 2(2X - 1) (3A + 1 - a) 



(90) 



Now we present an exphcit form of the function x entering the definition 
of fl2^ . and introduced in fl25]) . Due to Eqn f[T^ we have 



-AAx2 (A + 1/2)X^\ 
^ + AAx2 {A + 1/2) 



(92) 



Next, we give here an exphcit expression for the expression dx/ dx — 2x/ x. 
We compute 



ax _ 2x ^ AAx^ / 1 - a _ «AAx^ _ 

dx X X \ X yX 2yX ^ ' \ yX 



(93) 
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for X < Xrr, and 



^ _ 2x 
dx X 



x^Z 



\\x^ (\ - a 



Z \ X 
aXXx 



aY 



+ 



2yZ^ 



-2ycr + AAx^ {A +1/2) 



3 + ^|). (94) 



for X > Xcr- 

Note, that the function y, Eqn (1191) . is chosen so as to cancel the term 
— 2AAx/(A^(l + y"^)) arising from the last term in expression fl2^ (see the 
first term on the r.h.s. of flHHl) and the first term on the r.h.s. of f fT7|) ). With 
the help of Eqns fl29l) . fl93l) and fl88l) we obtain following expression for the 
function if) in the domain x < Xcr'- 



x^ 

^ = — 



iXXx A +1/2 



AX'^x 



2aXXx'^ 



X 



[i + y^yx 



+ 



2X^XXx^y {A - a) 

(1 + 7/2)2X2 



Y 
X 



XXx^ A + 1/2 



yX{l + y^ 

4(AA)Vt/ {A + 1/2) 



xY 



2x 



(1 + 2/ 



2X^ - XX 



2\2 



2 X 
' 4AV2(i _ y2 

X2 (1 + 2/2)2X2 



1 - Y A2AAx2 fax^ Y^ 

+ 77- 



X 



6ax 



2/2 X3 



2/1 + 



2?/ y 



3(1 



a] 



2(2Z- 1) (3A + 1 



X 



X2 



(95) 



Now, using Eqns (l92l) . ( l29l) . ( 193|) . (171), we find expression for function ■?/^ 
in the domain y > y^r- In fact, to obtain the equation on the parameter A 
we need to know only the part of in {?/ > ycr}-, proportional to A2. For 
this reason we write out only this part: 



,x^{2xY^ 2Y XXx' 



^ = 2X'-{—— + — + 



ax^ y2 go-j; 



y^\y Z"^ Z 2Z Xy"^ Z"^ 
2{2X - 1) (3A + 1 - « 



yZ 



3(1 - a) 



+ term proportional to AA. (96) 
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Finally, we show Eqn (l30ll which was stated in Section [3l Indeed, the 
definitions of Y and Z and the second equation in (fT9l) imply that for y > y^r 



Z 2a ' ^ 

where B := ln(AAa;^). Using these relations and the equation (!96l) . we arrive 
at the desired relation fl30l). 



Appendix 4 

In this appendix we compute the partial derivatives of energy E = E(a, P) 
w.r.to parameters a, (3. Using expression (l55ll . we obtain 

oo 


and 

oo 

dE 







(2) 


, 5y a 


(1) 


I 1/ ^1/ 


(1+ 7/2)2] 


[ at a/3 ' 




' ap a/? ' 







9„ 





\ dy d . 




1 5y a 


(1) 


: ^ y dy 


(1 + 2/2)2 1 


[dtda^ 




' ap aa ' 




' ' p"^ da 



2y l(ay\\(s^. ^» (33) 



Recall the notation 

X := IH and Z := 1 (99) 

2y 22/ ^ ^ 

and let 

X,^ := 1 + —-^ and Z,^ := 1 - ^— (100) 
From Eqn (fT9!) we find 

a?/ AAx^ 1 , , 
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and 



dy_ _ XX 1_ 

dye 



— X 



X > Xr 



XXxt 1 



dp 2(3 Xer 

Using Eqn fllOip . we obtain the time derivative of || in the leading approx- 
imation in 1/ In j • 



d f dy 
di \dp 



XXXx^ 



3 + 



aXXx^ Y 

~Wx 



X < Xr 



(102) 



d f dy 
di \M 



XXX [3 + 



2y2 Z 



2x1 



— X 



Qxt. 



Z V X(.:j^ J Xq^ 

In a similar way we find derivative of y w.r. to a: 

dy XXx^ In (|) 



X > Xcr- (103) 



da 



dy_ 
da 



XXl_ 



X 



xMn ( - ) + 2x' 

X 



X < Xr 



In 



Xr 



dye 



XXx"^ 



da 2 Xer 

Taking the time derivative of Eqn (11041) . we obtain 



d f dy 
dt V da 



XXXx^ ( XXx 



2AX2 



In 



AA 



(104) 



X X Q/f ^ 



y 



y^x^-^ 1+1/2 1- 



In 



X 



3 + 



oAAa^F 



X < Xer, (105) 
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d f dy 
dt \ da 



AAA J <J 2y z 



'2XZ{ 



2x' 



In ^ 



XXx^ In 



AA 



cr -lA 

+ 1/2 



-X In I - 



y\.cr 

X > Xcr- (106) 



Note that the main contribution to the partial derivatives 0971) and (1981) 
comes from the domain x ~ Xcr- Both derivatives are sums of terms propor- 
tional to A^ and to AA. The coefficients for these terms are of order of 1. As 
result, since we assumed that |AA| ^ A^, we have to find in the expressions 
for fl97j) and fl98l) only the terms proportional to A^. 

Using Eqns (110 II) - (11061) we can write the r.h.s. of ( 1971) and ( l98l) in a 
more explicit form 



IdE 



dxx^XX 



Y 



-X' 



dxx'^XX 



1 _|_ aXXx^ \ 
2y J 

Y 



3 + 



Y 



Ixp' \ _|_ aAAxf 



1x1 



X 



2?/ 

6x5?. 



X 



2x r 



2|/ 

«AAx^ 



2xy 

y 

Y 



2y^ 1 



a\\x'' 



+ 



1/ 1 



\\\x^ 
2?/ 



1 + 



— X 



22/c' 



(107) 



IDE 
1'da 



dxx^XXY ( XXx^ 




aXXx^ Y 



In 



/AA 



2y2 



lAAa;"' 



X 



6xL In ( 



1 + 



2?/ 



+ 



3 + 



2?/ 

aAAx^ Y 



+ ^ln f^)FKP 

?/ Vx 



M/2 
dxx^AAF 



2yM 1 - 



2j/ 



2y2 



2?/ 



2x1 In ^ 



1 + 



xMn(^ 

X 



X" 



1/(1-^) 



2?/cr - AAx^ In 



/AA 



1/2 



WM, September 15, 2009 



31 



2x Y /2x^rln(^|;; 




x'\n{^)\}. (108) 

-'cr 

Let a = a{P) be a solution of Eqn (I3U]) . We find numerically (see Ap- 
pendix 5) that P changes on the interval (0,/5o], where Po = 1.0405 (the 
corresponding value of a is ao = 0.65436.) Using expressions for |^ and 
derived above, we show numerically that the function 

is negative for (3 = (3q and for /3 0, with E{a, (3) having absolute minimum 
at /3 = [3q. 

Appendix 5 

Numerical calculations with help of Eqns ( ISTl) and ([5]) show that there is 
a point (tto, /^o), 

(3o = 1.0405, tto = .6543626, 

so that the equation (l50l) . I{a,P) = 0, has no solution for a < ag and for 
P > Pq. Moreover, the solution of the equation I{a,P) = for /3 determines 
a double- valued function P = /3(a), whose branches coalesce at a = cto and 
have different derivatives there (see Eqns flllOp and fillip below). Moreover, 
I{a,P) = has the unique solution Pq at a = ao. Hence the solution of the 
equation I{a,P) = for a also leads to a double- valued function a = a{P). 
Numerical calculations give the following expansions for the lowest branch, 

P = Po-pi{a-ao)-p2{a-aof, (110) 

a > ao, and and for the distance. A, between the branches along the a-axis, 

A = 7i(/5o-/5)-72(/3o -/?)', (Ill) 

where 

Pi = 2.54732, P2 = 13.8297, (112) 
71 = .08029, 72 = .42736. (113) 



(Solving fill op for a gives the lower branch of the function a = a{P). Adding 
fillip to this solution gives the upper branch of a = a{P).) 
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To find the second "end" point on the a— interval we check the point 
a = 1 where the dependence of ?/ on x in Eqn ( JT9l) can be found in an 
exphcit form. To do this we note that fl2T|) and fl22l) with a = 1 imply that 



7 



and XXyc 



(114) 



P ' " 2772 In 7 27 In 7 

We also have Ucr/xcr = 2/3. For a = 1 solvability condition of Eqn ( fT9l) is 



1/2 



Indeed, set 



7 > e 



y = T/cr^;, z = 1 + (5, 



X 



+ T 



In the range < 5 ^ 1 we have 

2 



1 



r-n^-2in7 



(115) 
(116) 

(117) 



From this equation we see that jS should satisfy the inequality given in Eqn 

(ma). 

Now we set y = UcrZ. For z < 1 we obtain from the first equation in (ITQll . 
with a = 1, and from flll4p the following cubic equation for the ratio 



X 
Vcr 



A f X \'^ln(72;) x 
— — -^-^ + z = Q. 

27 Vl/cr/ In 7 Vcr 



118) 



Solution of Eqn (11181) in the range z < 1 is 

1 \ 1/2 

m7 \ 



for 72; < 1 and 



X 

Vcr 



o , 1117 . 

— rsm 

m(72;) 



Vln(7^)/ln7 
3 yl — ln(72;)/ ln7 



(119) 



(120) 



for 72; > 1. 
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In the range z > 1 the ratio x/ycr solves the following cubic equation (see 
the second equation in ( fT9l) ) 

3 

+ 2-z = 0. (121) 



4 ln(7z) f X \ X 



cr 



27 In 7 \ycrj y, 
Let zq be solution of equation 

^ = ^,,.2/J^. (122) 

We split the semi-interval z > 1 into two sub-intervals. In the interval 
1 < z < Zq we have 

Vcr V ln(7z) 



where 



vr 1 , / ^l-(2-^)2 1n(7^)/ln7 \ ^ 



>^ H — arctan — = — , 1 < z < 2, 

6 3 V (2-;2)v/ln(7^)/ln7 J 

0=- + -arctan J v / ^ — \^2<z<zq. 124 

3 3 Vv/l-(^-2)2 1n(7z)/ln7; ^ ^ 

In the range z > z^we have 

- = |(Q^/^ + T^Q~^/^) (125) 

VcT 2 ln(7z) 

where 



Using Eqns flll9l) - fll26p we obtain with help of numerical calculations, that 
Eqn (150!) at a = 1 has solution only as (5 goes to zero. This means that a = 1 
is the second end point of the a— interval. 
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